A systematic study of the superdeformation of Pb isotopes with relativistic mean field 
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The microscopically constrained relativistic mean field theory is used to investigate the superde- 
formation for Pb isotopes. The calculations show that there exists a clear superdeformed minimum 
in the potential energy surfaces with four different interactions NL3, PKl, TMl and NLSH. The 
excitation energy, deformation and depth of well in the superdeformed minimum are comparable for 
four different interactions. Furthermore the trend for the change of the superdeformed excitation 
energy with neutron number is correctly reproduced. The calculated two-neutron separation energy 
in the ground state and superdeformed minimum together with their difi'erences are in agreement 
with the data available. The larger energy difference appearing in superdeformed minimum reflects 
a lower average level density at superdeformations for Pb isotopes. 
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Superdeformation(SD) of atomic nuclei is one of the 
most interesting topics of nuclear structure studies. Over 
the past two decades, many rotational bands associated 
with SD shapes have been observed in several regions of 
the nuclear chart 0] , with 85 SD bands observed in nuclei 
with 79 < Z < 84 (the A ~ 190 region) alone, where an 
impressive number of results has been obtained. Unfor- 
tunately, despite the rather large amount of experimen- 
tal information on SD bands, there are still a number 
of very interesting properties, which have not yet been 
measured. The characteristic examples are the spin, par- 
ity and excitation energy relative to the ground state of 
the SD bands. The difficulty lies with observing the very 
weak discrete transitions which link SD levels with lev- 
els of normal deformation (ND levels). Until now, only 
several SD bands have been identified to exist the transi- 
tions from SD levels to ND levels in the A 190 region: 
two bands in jjg j^^Q^ and one band in each of ^^^pb 
HH, i92pb 0, and is^Hg ^. Less precise measurements 
have been achieved in ^^^Hg and Pb ^9] following 
analysis of the quasi-continuum component of the decay. 
Recently, the measurement of the excitation energy of the 
yrast (lowest energy for a given spin) SD band in ^^Spb jg 
reported together with earlier measurements of the 
excitation energies of SD states in i94pb and ^^^Ph, al- 
lows a systematic study of the energy of the SD well in a 
single isotope chain. Many theoretical models have been 
employed to study these superdeformed states of atomic 
nuclei. The Strutinsky method with a Woods-Saxon po- 
tential the Hartree-Fock-Bogoliubov method with 
different mean fiel d p arameterizations [l^ IT^ Q| , and 
the cluster model [l5| have provided the predictions on 
the excitation energy of SD bands, where a gross trend 
of decreasing energy with decreasing neutron number is 
obtained for Pb isotopes, but the absolute energies as 
well as their differences are not consistently reproduced 
by these models as the analysis in Ref. Considering 
that the excitation energy and the well depth of the SD 



minimum are amongst the most important factors which 
affect the decay of the SD bands to the grou nd state, 
the relativistic mean field (RMF) theory [ij has also 
been recently applied to estimate the excitation energies 
and depths of well for SD bands. The excitation ener- 
gies and depths of well in the SD minima of i94jjg and 
i94pb have been predicted in considerable agreement as 
compared with experiment |l7l |. However, a systematic 
study of the energy of the SD well in a single isotope 
chain presented in I?,ef.|l(1| with the nonrelativistic theo- 
ries, has not been performed in a relativistic framework. 
In resent years, the RMF theory has gained consider- 
able success in describing many nuclear phenomena for 
the stable nuclei (18L 1191] as well as nuclei even far from 
stability It has been shown that the RMF theory 

can reproduce better the nuclear saturation properties 
(the Coester line) in nuclear matter |2ll |. present a new 
explanation for the neutron halo and predict a new 
phenomenon — giant neutron halos in heavy nuclei close 
to the neutron drip line 121^ , give naturally the spin-orbit 
potential, the origin of the pseudospin symmetry [3 113 
as a relativistic symmetry |26l H7, '28] and spin symme- 
try in the anti-nucleon spectrum ,29.J . and present good 
description for the magnetic rotation j^D], the collective 
multipole excitations l3l| , the identical bands in superde- 
formed nuclei j and the excitation energies relative to 
the ground of SD bands Hence, here we will report 
a systematic investigation of SD states for Pb isotopes in 
the microscopic quadruple constrained RMF theory with 
pairing treated by the BCS method, and show an excel- 
lent empirical manifestations of this SD structure in Pb 
isotopes including the evolution of the excitation energy, 
depth of well, deformation, and shell structure as well as 
the comparison with the ND states. 

The starting point of RMF theory is a standard 
Lagrangian density where nucleons are described as 
Dirac particles which interact via the exchange of var- 
ious mesons including the isoscalar-scalar a meson, the 
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isoscalar- vector uj meson and the isovector- vector p me- 
son. The effective Lagrangian density considered is writ- 
ten in the form: 



(1) 



where ip = '0^7*' and tp is the Dirac spinor. Other sym- 
bols have their usual meanings. 

The Dirac equation for the nucleons and the Klein- 
Gordon type equations for the mesons and the photon are 
given by the variational principle and can be solved by ex- 
panding the wavefunctions in terms of the eigenfunctions 
of a deformed axially symmetric harmonic-oscillator po- 
tential or a Woods-Saxon potential(3J|]. The details 
can be also found in Ref. ^'^d references therein. 

The potential energy curve can be calculated micro- 
scopically by the constrained RMF theory. The binding 
energy at certain deformation value is obtained by con- 
straining the quadruple moment {Q2) to a given value /i2 
in the expectation value of the Hamiltonian [35| , 



m = (H) 



1 



M2j 



(2) 



where is the constraint multiplier. 

For the nuclei studied in this paper, the deformed har- 
monic oscillator basis is taken into account and the con- 
vergence of the numerical calculation on the binding en- 
ergy and the deformation is very good. The converged 
deformations corresponding to different p,2 are not sen- 
sitive to the deformation parameter /3o of the harmonic 
oscillator basis in a reasonable range due to the large ba- 
sis. The different choices of /3o lead to different iteration 
numbers of the self-consistent calculation and different 
computational time. But physical quantities such as the 
binding energy and the deformation change very little. 
Thus the deformation parameter /3o of the harmonic os- 
cillator basis is chosen near the expected deformation to 
obtain high accuracy and low computation time cost. By 
varying /i2, the binding energy at different deformation 
can be obtained. The pairing is considered by the con- 
stant gap approximation (BCS) in which the pairing gap 
is taken as 12 /^/~A for even number nucleons. 

The calculated potential energy curves for i90-204p|-| 
are exhibited respectively in the figures d 12 E] and 
Hfor the interactions NLsjl^, PKljl^, TMlO and 
NLSHjs^, in which the energy of the ground state is 
taken as a reference. The and V are respectively the 
excitation energy and depth of well of the SD minimum 
as shown in the subfigure for ^^^Pb. Similar patterns are 
found for all the effective interactions. Most of the curves 
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FIG. 1: The potential energy curves for i90-204p^ obtained by 
the constrained RMF theory with the interactions NL3, where 
the Ex and V represent respectively for the excitation energy 
relative to the ground state of superdeformed minimum and 
the depth of well of superdeformed minimum. The ground 
state binding energy is taken as a reference. 



t PKl 



4 - 



> 
^ 16 




LU 



] iiii|i m i|i m: 





|iiii| Mn |ii 




I 1 1 1 1 1 1 I itI 1 1 1 I 1 1 1 1 1 I 1 1 itI I I 1 1 1 1 1 1 I 1 1 1 itI 




|llll| MM | 





-0.5 0.0 0.5 -0.5 0.0 0.5 -0.5 0.0 0.5 -0.5 0.0 0.5 
FIG. 2: The same as Fig.[Tl but with PKl 



display a clear SD minimum for the Pb isotopes, especial 
for the nuclei with more neutron number. For i90,i92p|-,^ 
the calculated SD minimum is not very obvious, even dis- 
appears for the interactions NLSH. For ^^°Pb, the RMF 
theory predicts a considerable high excitation energy rel- 
ative to the ground of SD bands and shallow well in the 
SD minimum in comparison with its neighboring nucleus 
^^^Pb, which implies it is difficult to come into being the 
stable SD state. For ^^^Pb, although the well is shal- 
low, the excitation energy is relatively lower, which in- 
dicates that it is relatively easy to form the SD state as 
observed in experiment. Start from the ^^"^Pb, the RMF 
theory predicts that the excitation energy increases with 
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FIG. 3: The same as fig.[T] but with TMl 
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FIG. 4: The same as Fig. d but with NLSH 



the increasing of the neutron number in company with 
the increasing of depth of well. So, the SD states can be 
still formed in these nuclei, which agrees the experimen- 
tal observations. However for the iV > 118 nuclei, as the 
excitation energy is too high to make it difficult to excite 
the SD states. It may explain why only the SD nuclear 
states between iV = 110 and iV = 116 are observed in 
the Pb isotope chain. Besides the success in describing 
SD states, the RMF theory predicts an interesting fea- 
ture in the ground state. The evolution of shape from the 
prolate to the oblate, and finally to the spherical shapes 
are found in the Pb isotope chain. From Figs. 1-4, the 
ground state of ^^°Pb is exhibited a coexistence of the 
prolate and the oblate with about 5 MeV stiff barrier 
against deformation. Begin from ^^"Pb, the ground state 
gradually moves toward the oblate side with smaller and 
smaller deformation with the increasing of neutron num- 



> 

CD 



LU 



1 


' 1 

-Q - - 


' 1 

NL3 


1 


1 


1 


I'M. 

p 

// p 




-0-- 


PK1 








— 


-A-- 


TMl 








p'// — 
,-;o' ,-' 




-V-- 


NLSH 










• 


EXP 




















-',o 


A' 










-;,o' _ 
1 






" 

-1 1 


1 1 


"■A-" 

1 1 


V.'.'- 
.-•A'' 
1 


1 


Pb J 
1 . 1 .- 


108 


110 112 


114 


116 


118 


120 122 



Neutron Number 

FIG. 5: The SD bandhead energies of Pb isotopes as 
the functions of neutron number obtained by the constrained 
RMF calculations with the interactions NL3 (open squares), 
PKl (open circles), TMl (up triangles) and NLSH (down tri- 
angles) in comparison with the experimental data (filled cir- 
cles) . 



ber. Finally well spherical ^"^Pb are seen. 

In Fig. 5, the calculated bandhead energies for SD 
bands as the functions of neutron number are displayed 
for Pb isotopes, where the open squares, open circles, 
up triangles, down triangles, and filled circles stand re- 
spectively for the RMF calculations with the interactions 
NL3, PKl, TMl, NLSH, and the data available [13 ■ 
First of all, trend for the change of the excitation en- 
ergies with neutron number are correctly reproduced for 
all the effective interactions. Especially for the NL3, the 
calculated excitation energies are in excellent agreement 
with the data. The maximum deviation between theory 
and experiment is less than 0.34MeV. For the TMl inter- 
actions, the excitation energy is overestimated for ^^°Pb 
and underestimated for i92,i94p|^ ^ jump appears in 
^^''Pb, which disagrees the experimental data. But the 
trend for the change of the excitation energies with neu- 
tron number is in accordance with the RMF predictions. 
For the PKl interactions, except the excitation energy of 
^^^Pb is overestimated, theoretical predictions is highly 
consistent with the experiment. For the NLSH interac- 
tions, the RMF theory fails to reproduce the SD mini- 
mum for i92,i94pj-|^ i^^-j. reproduce the SD minimum for 
^^^Pb. In addition, the trend for the change of excitation 
energies with the neutron number is consistent with that 
in the other interactions. Compared the RMF calcula- 
tions, the Strutinsky method with a Woods-Saxon poten- 
tial [llj, the Hartree-Fock-Bogoliubov method with dif- 
ferent mean field parameterizations [l^ , and the 
cluster model0 predict only a gross trend of decreasing 
energy with decreasing neutron number. The absolute 
energies and their differences are not consistently repro- 
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TABLE L The quadruple deformation /32 and the depth of 
the superdeformed minimum V in the superdeformed states 
of i90-204p|^ obtained by the constrained RMF theory with 
the interactions NL3, PKl, TMl and NLSH. 
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0.46439 


0.42181 


0.22936 






0.66217 


0.34837 


0.25777 




194p^^ 


1.37963 


0.29612 


0.59873 




196pb 


1.60176 


1.02856 


1.41779 


0.46892 


198pb 


2.17652 


1.66011 


1.56866 


1.51712 


200p^ 


2.83555 


2.34135 


1.75199 


2.17761 


202p^^ 


3.25852 


2.37082 


1.48711 


2.17321 


204p^ 


3.51464 


2.24769 


1.49029 


1.96709 



TABLE II: Two-neutron separation energy in the ground 
state and superdeformation minimum obtained by the con- 
strained RMF theory with the interactions NL3, PKl, TMl 
and NLSH, in comparison with data available 
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<52n,ND 


[Me V ] 
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iNL3 


FKl 


i Ml 


iN-LO-tl 


EAF 


no 


17.862 


16.915 


16.716 


17.984 


18.400 


112 


17.537 


17.207 


17.135 


17.444 


17.810 


114 


17.173 


16.885 


16.636 


17.255 


17.320 


116 


16.739 


16.440 


16.124 


16.899 


16.820 


118 


16.287 


16.050 


15.595 


16.652 


16.295 


120 


15.709 


15.241 


14.886 


15.150 


15.837 


122 


16.360 


15.287 


14.927 


15.611 


15.318 






'S'2n,SD 


[MeV] 






N 


NL3 


PKl 


TMl 


NLSH 




110 


18.364 


18.025 


17.908 






112 


17.031 


17.334 


18.436 




17.16(4) 


114 


15.968 


16.266 


16.037 




16.31(4) 


116 


15.287 


15.081 


14.820 


15.499 




118 


14.648 


14.303 


14.051 


14.653 




120 


13.574 


13.161 


13.202 


13.361 




122 


12.906 


12.484 


12.601 


12.637 





duced by these model. It shows that the RMF theory 
gives a better description of the SD excitation energies 
for the Pb isotopes. 

Besides the excitation energy, the deformation and well 
depth of SD minimum are another two important param- 
eters which reflect the properties of the superdeformed 
states. In particular, the well depth affects the life time 
of the superdeformed states. In Table the quadru- 
ple deformation and the depth of the superdeformed 
minimum V in the superdeformed states are listed re- 
spectively in the upper and lower panels for i90-204pj^ 
Except for several exceptions, the deformation in SD min- 
ima lies systemically between 0.5 and 0.7 for four differ- 
ent interactions NL3, PKl, TMl and NLSH, which agrees 
the observation of superdeformed nuclei for excited states 
adopting ellipsoidal shapes with an axis ratio around 2:1 
[y. The RMF theory predicts for the height of the bar- 
rier is lower than 1 MeV for i90:i92p|-,^ higher than 1 
MeV for i96-2i4p^^ ^^le interactions. For ^^''Pb, the 
estimated barriers are considerable different from differ- 
ent interactions. 

Two- neutron separation energy defined as S2n{Z, N) = 
E{Z, N) — E{Z, N — 2) is sensitive quantity to test a mi- 
croscopic theory, where E{Z, N) is the binding energy of 
nucleus with proton number Z and neutron number iV. 
In Table m the two-neutron separation energies in the 
ground state 5*2™, nd and SD minimum S'2n,SD are shown 
respectively in upper and lower panels for the Pb iso- 
topes in comparison with the datalTnlliot. From there, 
it is found that the RMF calculations with four differ- 
ent interactions well reproduce the experimental data for 
'S'2n,ND- The maximum deviation between the calcula- 



TABLE HI: Two-neutron separation energy difference in the 
ground state and superdeformation minimum obtained by the 
constrained RMF theory with the interactions NL3, PKl, 
TMl and NLSH, in comparison with data available 







AS'2n,ND 


[MeV] 






N 


NL3 


PKl 


TMl 


NLSH 


EXP 


110 


0.325 


-0.292 


-0.419 


0.540 


0.590 


112 


0.364 


0.322 


0.499 


0.189 


0.490 


114 


0.434 


0.445 


0.512 


0.356 


0.500 


116 


0.452 


0.390 


0.529 


0.247 


0.525 


118 


0.578 


0.809 


0.709 


1.502 


0.458 


120 


-0.651 


-0.046 


-0.041 


-0.461 


0.519 






AS'2n,SD 


[MeV] 






N 


NL3 


PKl 


TMl 


NLSH 


EXP 


110 


1.333 


0.691 


-0.528 






112 


1.063 


1.068 


2.399 




0.85(8) 


114 


0.681 


1.185 


1.217 






116 


0.639 


0.778 


0.769 


0.846 




118 


1.074 


1.142 


0.849 


1.292 




120 


0.668 


0.677 


0.601 


0.724 





tions and data is less than 1.7 MeV, especially for the 
NL3, the deviations is within 1 MeV. For S2n,su, the 
calculations with four different interactions are compara- 
ble and close to the data available. Both the calculations 
and experiment show that the 5'2n,ND and S'2n,ND vary 
smoothly with the neutron number. No sharp drop in 
the binding energy is seen from the S2ni which indicates 
no significant shell gap appearing in the Pb isotope chain 
whether the ground state or SD states. 

In order to reveal further the detailed information on 
shell structure, the two-neutron separation energy differ- 
ences AS2n{Z, N) = S{Z, N)-S{Z,N + 2) are presented 
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in Table UlTI where the AS'2„,nd and AS2n,SD represent 
respectively for those in the ground state and superde- 
formed minimum with the data for comparison ^3 E3- 
From the Table IIIII it is found that the experimental 
AS'2ti,nd changes typically around O.SMeV, is well repro- 
duced in the RMF calculations for all the interactions 
with several exceptions. Compared with the other inter- 
actions, the NL3 gives better agreement with experiment. 
Only the deviation is relatively large at ^°^Pb. Further- 
more, the calculated AS'2n,ND shows very little differ- 
ences for these nuclei with neutron number from A'^ = 112 
to = 118, is consistent with experiment data, which 
suggests without a shell closure appearing at TV = 112 or 

= 114 predicted in other calculations. Compared with 
the ND states, the SD separation energies S2n are sig- 
nificantly larger than the typical ND value of 0.5 MeV, 
possibly reflecting a lower average level density at su- 
perdeformations. In particular, the AS'2n,sD presents a 
obvious difference for different nuclei. The A5'2ra,SD for 
these nuclei with N = 112, 114, 118 is much larger than 
that for their neighboring nuclei, suggest a larger shell 
gap in these SD states of "a.iw.nspi^ 

In summary, the superdeformation in i90-204p|-j jg 
vestigated by the microscopic quadruple constrained rel- 
ativistic mean field theory with all the most used interac- 
tions, i.e., NL3, PKl, TMl and NLSH. The calculations 
show a clear SD minimum at nearly all the potential en- 
ergy curves for the Pb isotopes with similar patterns for 
all the effective interactions. Trend for the change of the 
excitation energies with neutron number are correctly re- 
produced. The calculated deformation in SD minima lies 
systemically between 0.5 and 0.7, is consistent with the 
observation of experiment. The two- neutron separation 
energies in the ground state and the SD minimum are 
well reproduced with varying smoothly with the neutron 
number. Compared with the ND states, the SD separa- 
tion energies 6*2,1 are significantly larger than the typical 
ND value of 0.5 MeV, possibly reflecting a lower average 
level density at superdeformations. 
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Foundation in University of Anhui Province in China. 



* Electronic address: 'jianyou@ahu.edu.cnl 
[1] B. Singh, R. Zywina, and R. B. Firestone, Nuclear Data 

Sheets 97, 241 (2002). 
[2] T. L. Khoo, M. P. Carpenter, T. Lauritsen et al., Phys. 

Rev. Lett. 76, 1583 (1996). 
[3] G. Hackman, T. L. Khoo, M. P. Carpenter et al., Phys. 

Rev. Lett. 79, 4100 (1997). 



[4] A. Lopez-Martens, F. Hannachi, A. Korichi et al., Phys. 

Lett. B 380, 18 (1996). 
[5] K. Hauschild, L. A. Bernstein, J. A. Becker et al., Phys. 

Rev. C 55, 2819 (1997). 
[6] A. N. Wilson, G. D. Dracoulis, A. P. Byrne et al., Phys. 

Rev. Lett. 90, 142501 (2003). 
[7] S. Siem, P. Reiter, T. L. Khoo et al, Phys. Rev. C 70, 

014303 (2004). 

[8] T. Lauritsen, T. L. Khoo, I. Ahmad et al., Phys. Rev. C 

62, 044316 (2000). 
[9] M. S. Johnson, J. A. Cizewski, K. Y. Ding et al., Phys. 

Rev. C 71, 044310 (2005). 
[10] A. N. Wilson, A. K. Singh, H. Hiibel, et al., Phys. Rev. 

Lett. 95, 182501 (2005). 
[11] W. Satula et al., Nucl. Phys. A529, 289 (1991). 
[12] S. J. Krieger et al., Nucl. Phys. A542, 43 (1992). 
[13] P.-H. Heenen, J. Dobaczewski, W. Nazarewicz et al., 

Phys. Rev. C 57, 1719 (1998). 
[14] J. Libert, M. Girod, and J.-P. Delaroche, Phys. Rev. C 

60, 054301 (1999). 
[15] G. G. Adamian, N. V. Antonenko, R. V. Jolos et al., 

Phys. Rev. C 69, 054310 (2004). 
[16] B. Serot and J.D. Walecka, Adv. Nucl. Phys. 16, 1 (1986). 
[17] G.A. Lalazissis, P. Ring, Phys. Lett. B 427, 225 (1998). 
[18] P. -G. Reinhard, Rep. Prog. Phys. 52, 439(1989). 
[19] P. Ring, Prog. Part. Nucl. Phys. 37, 193(1996). 
[20] J. Meng, Nucl. Phys. A 635, 3(1998). 
[21] R. Brockmann and R. Machleidt, Phys. Rev. C 42, 

1965(1990). 

[22] J. Meng and P. Ring, Phys. Rev. Lett. 77, 3963(1996). 
[23] J. Meng and P. Ring, Phys. Rev. Lett. 80, 460(1998). 
[24] A. Arima, M. Harvey, and K. Shimizu, Phys. Lett. B 30, 
517(1969). 

[25] K. T. Hecht and A. Adler, Nucl. Phys. A 137, 129(1969). 

[26] J. N. Ginocchio, Phys. Rev. Lett. 78, 436(1997). 

[27] J. Meng, K. Sugawara-Tanabe, S. Yamaji, and A. Arima, 

Phys. Rev. C 59, 154 (1999). 
[28] J.-Y. Guo, R.-D. Wang, and X.-Z. Fang, Phys. Rev. C 

72, 054319 (2005) 
[29] S. G. Zhou, J. Meng, and P. Ring, Phys. Rev. Lett. 91, 

262501 (2003). 

[30] H. Madokoro, J. Meng, M. Matsuzaki, S.Yamaji, 

Phys.Rev. C 62, 061301(2000). 
[31] Z. Y. Ma, A. Wandelt, N. V. Giai, D. Vretenar, P. Ring, 

L. G. Cao, Nucl. Phys. A 703 222(2002). 
[32] J. Konig and P. Ring, Phys. Rev. Lett. 71, 3079(1993). 
[33] Y. Gambhir, P. Ring, A. Thimet, Ann. Phys. (N.Y.) 198, 

132 (1990). 

[34] S.G.Zhou, J.Meng, P.Ring, Phys. Rev. C 68 (2003) 
034323 

[35] P. Ring and P. Schuck, The Nuclear Many Body Problem, 

(Springer 1980). 
[36] G. A. Lalazissis, J. Konig, and P. Ring, Phys. Rev. C55, 

540 (1997). 

[37] W. H. Long, J. Meng, N. V. Giai, and S. G. Zhou, Phys. 

Rev. C 69, 034319 (2004). 
[38] Y. Sugahara, and H. Toki, Nucl. Phys. A 579, 557 (1994). 
[39] M.M. Sharma, M.A. Nagarajan, and P. Ring, Phys. Lett. 

B312,377 (1993). 
[40] G. Audi and A. H. Wapstra, Nucl. Phys. A 595, 

409(1995). 



